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Abstract 



This paper presents the vortical and self-similar solutions for 2D compressible Euler equa- 
tions using the separation method. These solutions complement Makino's solutions in radial 
symmetry without rotation. The rotational solutions provide new information that furthers 
our understanding of ocean vortices and reference examples for numerical methods. In ad- 
dition, the corresponding blowup, time-periodic or global existence conditions are classified 
through an analysis of the new Emden equation. A conjecture regarding rotational solutions 
ON ' in 3D is also made. 
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c3 ! 1 Introduction 

In fluid mechanics, the N-dimensional isentropic compressible Euler and Navier-Stokes equations 
are expressed as follows: 

p t +V ■ (pu) =0 

(1) 

p{u t + (u ■ V)u\ + K\7p~< = 0, 
where p — p(t,x) denotes the density of the fluid, u = u(t,x) = (u\,U2, ■ ■ ■ , uat) S R n is the 
velocity, x = (x±, x-i, ■ ■ ■ , xn) £ R N is the space variable, and K > 0, 7 > 1 are constants. 
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The Euler equations have been studied in great detail by numerous scholars because of their 
significance in a variety of physical fields, such as fluids, plasmas, condensed matter, astrophysics, 
oceanography, and atmospheric dynamics. These equations are also important in physics and 
are widely used in different areas of study. For instance, the Euler system is the basic model 
for shallow water flows [3]. It also provides a good model of the superfluids produced by the 
Bose-Einstein condensates in the dilute gases of alkali metal atoms, in which identical gases do 
not interact at very low temperatures [3]. At the microscopic level, fluids or gases are formed by 
many tiny discrete molecules or particles that collide with one another. As the cost of directly 
calculating the particle-to-particle or molecule-to-molecule evolution of fluids on a large scale is 
expensive, approximation methods are needed to considerably simplify the process. Therefore, at 
the macroscopic scale, the continuum assumption, which considers fluids as continuous, is applied 
to the modeling. Here, the Euler equations provide a good description of the fluids at the statistical 
limit of a large number of small ideal molecules or particles by ignoring the less influential effects, 
such as the self-gravitational forces and relativistic effect [T]. For a mathematical introduction of 
the Euler equations, readers are referred to [5] and [2J. 

The construction of analytical or exact solutions is an important area in mathematical physics 
and applied mathematics, as it can further classify their nonlinear phenomena. For non-rotational 
flows, Makino first obtained the radial symmetry solutions for the Euler equations (fTJ in R N in 1993 
[5] . A number of special solutions for these equations [5] , [7] , [T2J , [13] , and [T3] were subsequently 
obtained. For rotational flows, Zhang and Zheng [TS] constructed explicitly rotational solutions 
for the Euler equations with 7 = 2 in 1997: 

P= s^' = u * = h ix - y) > (2) 

where x = r cos 9 and y = r sin 9. 

Very recently, Kwong and Yuen [5] constructed a family of rotational solutions for the Euler- 
Poisson equations 

{p t +V ■ (pu) =0 
p{u t + (u ■ V)u)+KVp =-pV$ (3) 
A$(t,f) = 27rp, 

with N = 2 and 7 = 1: 

= H&re'C'M')), u{t,x) =^ {x ,y)+^(-y, x ) 
fi (*) = + = «o > 0, 6(0) = 01 (4) 

e /(S) = #, HO) = a, /(0) = 0, 
with arbitrary constants £ ^ 0, ao, a\ and a. 

The rotational case complements Yuen's solutions without rotation (£ = 0) [TT]. In this paper, 
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based on the foregoing development, we provide the corresponding vortical flows for 2D compress- 
ible Euler equations ((TJ with 7 > 1 in the following result. 

Theorem 1 For 7 > 1, there exists a family of vortical flows in radial symmetry for the com- 
pressible Euler equations (QP in 2D, 

»«((-^«+«)^' ) 

f ) = — 

u(t,x)=^(x,y)+-^(-y,x) ^ 
k a (*) ^W^t = 0, o(0) = a > 0, 6(0) = ai, 

2. 2 

wzi/i the self-similar variable s = x a ^2 and arbitrary constants £ ^ 0, ao, ax, and a. 

Remark 2 This result complements Makino's solutions in radial symmetry without rotation = 
0). The vortical solutions provide new information that may further our understanding of 
oceans vortices and reference examples for numerical methods in computational physics. 

Remark 3 Zhang and Zheng's solution ^ for 7 = 2 is a special case in our solutions {5p. 



2 2D Vortical and Self-similar Flows 



Here, we provide a lemma for the vortical flows in 2D for the mass equation (HJi. This lemma 
originates in Kwong and Yuen's paper [5], which constructs periodic and spiral solutions for 2D 
Euler- Poisson equations §S§ with 7 = 1. 



Lemma 4 ([5]) For the equation of the conservation of mass, 

p t +V- (pu) = 0, 

there exist the following solutions: 

p(t, x) = p(t, r) = I ' , u (, x) = — (x, y) + -y, x) 

a(t) z a(t) r 



with the radial r = \J x 1 + y 2 and arbitrary functions f > 0; G and a(t) > 6 C 1 . 
Proof. The functional structure 

p(t,x) = p(t,r), u(t,x) = (x,y) H (~V,x) 

r r 



(6) 



(7) 



(8) 



with an arbitrary C 1 function F(t, r), can be plugged into the mass equation to verify the result: 

(9) 
(10) 





Pt + v • 


(pu) 




d 
dx 


Fx Gy' 
p p — 

r r 


d 
dy 


Fy Gx' 
p hp 

r r 
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d Fx f d Fx\ d Gy . d Gy . 

pt + inP— +p (in— - (^zP)— - Pirn—) 



dx r \dx r J dx r dx r 
d Fy f d Fy\ d > Gx . d Gx 

Yy 9 — + P [aUTJ + { d-y P) — + ^ V } (U) 



X £ X ( X\ X f X 

Pt + Pr h P F r - - + p pFx— 

r r \ r J r r r° 

xGy xy x y Fy { „ V\ V 

Pr pG r h pGy^r + p r h P(F r ~ - 

/r* rr* rv* <v* tv* tv* \^ f* / f* 



F y y Gx / y\ x y 
p pFy— + p r hp G r - pGx— (12) 



x Fx / x\ x F x 

Pt + Pr 1- P F r - - + p pFx- 

r r \ r / r r i 

y Fy . ( „ y\y . f y 

Pr +p[F r -)-+p pFy-r (13) 

r r Vr/r r r a 

1 , * 

= p t + p r F + pF r + pF-. 14 
r 



Then, the self-similar structure is taken for the density function, 



p(t,x)=p{t,r) = f -^^ (15) 

id F(t, r) = ^pyf for velocity u to balance equation \1 If : 

d Aggj) d f(-^Tj) a(t)r AogjJaft) AggOaft) , . 

~ di a(t) 2 ih a{tf a(t) a(t) 2 'ajtj a(t) 2 a(t) { ' 

-2a(t)f(r/a(t)) a(t)r f (r / a(t)) 



a(t) 3 a{ty 
f(r/a(t))a(t)r , f(r/a(t))a(t) , f(r/a(t))a(t) 



(17) 



a(t) 3 a(t) a{t) 2 a(t) a{t) 2 a(t) 

= 0. (18) 

The proof is completed. ■ 

The computational proof for Theorem 1 is as follows. 
Proof of Theorem 1. The procedure for the proof for vortical fluids is similar to that for 
non- vortical fluids [5] , [2] . It is clear that the following function 

p(t,x) = 4S> =^( x >y)+-!L(-y> x )> ( 19 ) 

a(ty a(t) a (t) 

with an arbitrary C 1 function / > of the self-similar variable s = x a t\yi = a ( t yt an d a unde- 
termined C 2 time-function a(t) > 0, satisfies Lemma |4] for the mass equation ([T])i . For the first 
momentum equation (jT])2i , we obtain: 

f d du\ dui\ , , 3 

= p {dt^ + u ^ + u ^) +K ^ Tx p (20) 
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1 du\ 



U2 



dui 
dy 



f(s) 



7-1 



1 d 



a(t) 2 (7-i) a(t) 2 dx 
f{sV~ 2 2x 



f(s) 



dt Ul + Ul ~te +U *W + K \(t)^-V a {tf 



at yait)^ 



w?y 



d_ ( a(t) , 



= P 



r ait)*) dy \a(t) X a 

I d(t) a 
\ a ( t ) a 



\wj x ~s^y) dx \a(t) x 7rfty?y 
i 



■v 



/«2 



off) 
a(t) 



«W 

<z(i) a(ty 



x + if 7 



, 7 -2 



px 

^tj 2 ^ 



A + 2if 7 /( S )^ 2 /(s) 



»(t)*r 

7-2 



■2x/(s) 



with the Emden equation 



a(0) = ag > 0, d(0) = a\, 



with arbitrary constants £ and A > 1. 

To ensure that equation (|26p is zero, the following ordinary differential equation is required: 

A + 2K 1 f(s)'<- 2 f( S ) = 
/(0) = a > 0, /(0) = 0. 

The function /(s) in the foregoing ordinary differential equation can be solved exactly: 

A( 7 -l) , 



2K-) 



For the second momentum equation (PJ22, we also have 

= P 
d 



= P 



d du 2 du 2 \ „_ x d 

-7TU 2 + Ui— h U 2 — — + A7// 7 — p 

at ox ay ) ay 

2y 



du 2 

—u 2 + Ui—— 

ot ox 



u 2 



du 2 
dy 



7-2 



ait) 2 ^- 1 ) a(ty 



dt V^(FF + a{t)y ) + \a(t) 

my] JL ( _£_ 

a{t)V) dy \a(ty 
+ \a(t) a(t) 2 J y + \a 



2go(t) , 
a (t)3 ■ 



a{t)V 

UI 
(*)' 



+ Kl^^2yf{s) 



i 

«(() 2 



/to" 1 



a(t) 



\a(i) a(ty 



y + K,±U^2yf(s) 



py 

a(t) 2 ^ 



A + 2K 1 f{ S y- 2 f{ S ) 
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= 0. (36) 

with the Emden equation (f2"T|) and function (|29p. 

Therefore, for £ ^ 0, we have the vortical and self-similar flows for the compressible Euler 
equations {J) in 2D [9]. 

To ensure the non-negative density function, we require that 

max ( (- Xi -2K^ s + a ) ^ 1 ' 
P(*,50 = — ^ ^ ^ (37) 

This completes the proof. ■ 
Lemma 5 For the Emden equation, 



a(0) = ao > 0, d(0) = oi, 
wiift arbitrary constants £ ^ 0, A and 7 > 1, the total energy is 



(38) 



- — +2^)2 + (27 _ 2)a(t)27 - 2 - ( 39 ) 

VFe /iaue the following. 

(1) For 1 < 7 < 2, if E(0) < 0, the solution is time-periodic; otherwise, the solution is global. 

(2) For 7 = 2, 

(2a) with £ 2 > —A, the solution is global; 

(2b) with £ 2 < —A, the solution blows up at a finite time if 



Ox < ^- X ~ e - (40) 
ao 

otherwise, the solution is global. 
(3) For 7 > 2, 

(3a) with A > 0, the solution is global; 
(3b) with A < and a constant a,Max — 
(3bl) and a > a M ax, 

if E(0) < F pot (aMax) or E(0) > F pot {aMax) with a\ > 0, the solution is global; otherwise, the 
solution blows up at a finite time. 
(3bII) and a < a Max , 

if E(0) > F pot {aMax) with a\ > 0, the solution is global; otherwise, the solution blows up at a finite 
time. 

Proof. For equation f|38[) . we multiply a(t) and then integrate it, as follows: 

— + 2aW + (2 1 -2)a(t)^- E{ ^ (41) 
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with a constant E(Q) = ^ + 



A 



2 1 2a 2 1 "•- o\ 



2af t (2 7 -2)a^ 

We define the kinetic energy as: 



a{tf 



and the potential energy as: 



F kin ■= (42) 



Fpot ~ 2a(t) 2 + (2 7 -2)a(t) 2 7-2- (43) 
The total energy is conserved thusly: 

dE(t) d , 

-^ = Tt (F kin+ F pot ) = 0. (44) 

The classical energy method for second-order autonomous ordinary differential equations (which 
readers can refer to pages 793-798 in [ID]), can be applied to analyze the corresponding qualitative 
properties of the Emden equation (l38l) . 

(1) For 1 < 7 < 2, there exists a unique global minimum for the potential function (|43[) for a(t) > 0, 
and lim F pot (a(t)) — +oo and lim F pot (a(t)) = 0. If 

a(t)-s-0+ a(t)-s- + oo 

E <°>=f + lr (2T - 2 V' <0 - w 

the solution is time-periodic; otherwise, it is global. 

(2) For 7 = 2, the ordinary differential equation (|3"8"f degenerates into the classical Emden equation 



& W = ^F- (46) 
a{t) 3 



(2a) For ^ 2 > —A, the solution is global. 

(2b) For £ 2 < —A, the potential function (|4"3"]l is strictly increasing with lim F pot {a{t)) = 0. 

a(t)— »+oo 

The solution blows up at a finite time, if 



a, < =-3; (47) 
a 

otherwise, the solution is global. 
(3) For 7 > 2, 

(3a) with A > 0, the potential function (1431) is a decreasing function of a(t) > 0. Thus, the solution 
is global. 

(3b) With A < 0, the potential function (|43[) achieves a unique global maximum, F pot (aMax), with 
a-Max = (^t) 2 '' 4 for a(t) > 0. 
(3bl) And a Q > a Max , 

if E(0) < F pot {aMax) or -E'(O) > F pot (aMax) with a\ > 0, the solution is global; otherwise, the 
solution blows up at a finite time. 
(3bII) and a < a Max , 
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if E(0) > F pot (aMax) with a\ > 0, the solution is global; otherwise the solution blows up at a finite 
time. 

The proof is now complete. ■ 

The foregoing lemma makes it easy to determine the blowup or global existence of the corre- 
sponding solutions J5]) for the compressible Euler equations JI]) in 2D. 

Corollary 6 The total energy is defined 

E(t) - — + 2^)2 + (2 7 - 2)a(i)27-2 > ( 48 ) 

for the Emden equation 

(1) For 1 < 7 < 2, if E(0) < 0, solution {3P is time-periodic; otherwise, solution |5|] is global. 

(2) For 7 = 2, 

(2al) with £ 2 > -A or £ 2 = -A and eti > 7 solution {5[) is global; 
(2aII) with £ 2 = —A and a\ < 0, solution [5$ blows up at T — —ai/ag. 
(2b) with £ 2 < —A, solution |2P blows up at a finite time if 



ai < VE0! ; (49) 
a 

otherwise, solution (0) is global. 
(3) For 7 > 2, 

(3a) with A > 0, solution ^ is global; 

i 

(3b) with A < and a constant aMax = ( -jr J > 
and a > aMax, 

if E(Q) < FpotiflMax) or E(0) > F pot {aMax) with a\ > 0, solution f5|) is global; otherwise, solution 
(0) blows up at a finite time. 
(3bII) and a < a Max , 

if E(0) > F pot (aMax) with a\ > 0, solution $5$) is global; otherwise, solution |3P blows up at a finite 
time. 



3 Conclusion and Discussion 

This paper provides a class of self-similar vortical flows for the 2D compressible Euler equations. 
The result presented herein complements Makino's solutions in radial symmetry without rotation. 
In addition, the corresponding blowup or global existence conditions arc classified by analyzing 
the new Emden equation Solution ([5]) is also the solution of the compressible Navier-Stokes 

equations in 2D: 

Pt+y ■ (pu) =0 

(50) 

p[u t + (u ■ V)n] + KVp' 1 = pAu, 
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with a positive constant fj,. 

Based on the existence of the rotational and self-similar solutions ([5]) in 2D, it is natural to 
conjecture that there exists a class of rotational solutions for the compressible Euler equations ([1]) 
in 3D, that complements those in radial symmetry, [5], |14j : 



(51) 



where 




ail n a k 

\k=i 



o»(Q) = am > 0, d,(0) = cii 



(52) 



with 

s = ^(^i)2. (53) 

fc=i flfc 

Here, £,d* ,d* 1; a > are constants, and functions ai — Oj(t).' 

~ t^tt /or i= 1,2,3 



(54) 



arbitrary constants a.;o and a.;i- 
Further research will be carried out to investigate the possibility of these vortical solutions in 3D. 
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